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Ionic Cloud Distribution close to a Charged Surface in the Pres- 
ence of Salt 
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Abstract. - Despite its importance, tlie understanding of ionic cloud distribution close to a charged macroion 
under physiological salt conditions has remained very limited especially for strongly coupled systems with, 
for instance, multivalent counterions. Here we present a formalism that predicts both counterion and colon 
distributions in the vicinity of a charged macroion for an arbitrary amount of added salt and in both limits of 
mean field and strong coupling. The distribution functions are calculated explicitly for ions next to an infinite 
planar charged wall. We present a schematic phase diagram identifying different physical regimes in terms of 
electrostatic coupling parameter and bulk salt concentration. 



Introduction. - Electrostatic interactions play a key role 
in controlling solubility, structure and phase behavior of 
macroions in aqueous solutions [1-3]. Examples of biolog- 
ically relevant macroion systems are charged lipid bilayers 
such as those found in mitochondrial membranes that contain 
considerable amounts of anionic cardiolipins or plasma mem- 
branes rich in anionic phospholipids {e.g. phosphatidylserines), 
stiff {e.g. DNA) or flexible {e.g. RNA) polyelectrolytes con- 
taining dissociated negatively charged phosphate groups, or 
charged polypeptides with a net charge depending on the disso- 
ciation equilibrium of various peptide moieties, as well as their 
complexes as encountered in the context of gene therapy [4] or 
self-assembly of viruses [5]. Instead of residing exactly on the 
charged surface, in order to minimize their electrostatic inter- 
action energy, the counterions needed to neutralize these sys- 
tems are distributed some distance away as a consequence of 
their translational entropy as well as the screening effects due 
to residual co- and counterions. It is the nature of the spatial 
distributions of all these various mobile ion species in the vicin- 
ity of a charged macroion that presents the biggest challenge in 
understanding charged (bio)colloidal systems [1]. This is the 



problem that we scrutinize in what follows. 

The traditional approach to charged macromolecular systems 
under salt-free conditions has been the Poisson-Boltzmann 
(PB) formalism, in which the Coulombic interaction between 
counterions is handled on a mean-field level [2,3]. However, 
in many biologically relevant situations the PB approximation 
breaks down; examples include most prominently multivalent 
counterions and highly charged surfaces. The most dramatic 
indication of this breakdown is the existence of attractive in- 
teractions between like-charged macroions [6-9] which on the 
PB level are known to be repulsive [10]. Consequently, there 
have been a number of attempts to assess corrections to the PB 
theory using, e.g., correlated density fluctuations around the 
mean-field distribution or additional non-electrostatic interac- 
tions [6, 11-14]. An alternative approach has been pioneered 
by Rouzina and Bloomfield [15] and elaborated by Shklovskii 
et al. [16], and later by Netz et al. [1, 17-19]. This approach 
leads to a new description of a system composed of a charged 
macroion and mobile counterions called the strong coupling 
(SC) theory. This description was shown to become exact in 
the limit of high surface charge, multivalent counterions, or 
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low temperature [17], clearly opposite to the PB limit, asymp- 
totically valid in the limit of low surface charge, monovalent 
counterions, or high temperatures. The PB and the SC theories 
thus asymptotically embrace all possible scenarios in the no- 
salt case. The SC theory has been applied with notable success 
to the case of charged macroscopic surfaces of various geome- 
tries with counterions [19]. 

However, under physiologically relevant conditions, the situ- 
ation is considerably more complicated. Biological systems al- 
ways contain significant concentrations of excess salts, which 
affect or quite often even govern their behavior. Overall, in 
biological conditions the reservoir salt concentration is typi- 
cally of the order of [Na^] = 100 mM. It is thus obvious 
that the average separation between nearest salt ions is small, 
implying that there is no justification to disregard the effects 
of salt, as opposed to the effects of counterions, in a consis- 
tent statistical-mechanical treatment of such systems. The tra- 
ditional approach in this context has been the Debye-Hiickel 
(DH) theory, related to the PB formalism, that applies in the 
hmit of small overall charges in the system [3, 11,20]. Apart 
from this limiting case, a coherent theoretical description in- 
cluding the SC limit has been missing. 

In this work we propose a consistent and systematic ap- 
proach to charged systems, composed of fixed macroions, mul- 
tivalent counterions (of charge valency +qc) as well as added 
salt (of cationic and anionic charge species of valency +g+ and 
— q_) in chemical equilibrium with a bulk reservoir. We con- 
sider the influence of added salt on ion distributions on the PB 
level as well as in the SC limit. The results are shown to fill 
the gap between the no-salt and high salt limits. In particular, 
we show how certain divergencies may be removed and nor- 
malizable ion density profiles may be obtained in the SC limit. 
Explicit calculations are carried out for ions at an infinite planar 
charged wall, where we also present subleading corrections to 
the SC results. 

Before we formulate the partition function for an interacting 
system of this type in a general form, let us introduce the rel- 
evant length scales and parameters of the problem, and show 
by scaling arguments what one should expect from a more rig- 
orous theory. First, we focus on counterions. An important 
length scale in the problem is the Gouy-Chapman (GC) length 
^ = \/{2Tiqc(JslB)- It measures the typical distance from a 
charged macroion surface (of surface charge density — (Ts) at 
which the electrostatic potential energy of a counterion inter- 
acting with the surface matches the thermal energy k-gT. Here, 
/b = / [Anek-QT] is the Bjerrum length, the distance at which 
the interaction between two unit charges equals thermal energy; 
in water I-q — 0.7 nm. The ratio between these two length 
scales yields an important dimensionless parameter, namely, 
the electrostatic coupling parameter [17] 
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embodying the relative contribution of ion-ion vs. ion-surface 
interactions. A related length scale is the lateral distance 
between counterions forming a strongly correlated quasi-2D 
layer close to a charged surface in the SC regime [1, 19], i.e. 



a± = \J [qc/T^Os) ~ V2Ei^i. This lateral distance is important, 
because in the SC regime we can express the concentration of 
the counterion layer as ric = as/{qcP-) = VI^'^Im) [15]. The 
counterion concentration within this strongly correlated diffuse 
layer should be independent of the bulk ion concentration, as 
long as the ionic strength, defined as / = i J2i=± ^ilh si§" 
nificantly smaller than Uc, i.e. I ^ Uc- Here and qi are the 
fugacities (to be defined later) and valencies of reservoir salt 
ions. This means that the counterion concentration on a surface 
should be much larger than the reservoir salt concentration. 

In a system that contains added salt there exists of course 
another independent length scale, the Debye screening length 
'dh — 1/ (Stt/Zb ) being related to the ionic strength /. It mea- 
sures the distance at which the Coulombic interaction between 
two point charges is screened out [20]. The effects of salt thus 
become unimportant when np, <C 1, where k = l^]^ denotes 
the inverse Debye screening length. Later on we will see that 
it is really the coupling constants nji and S that uniquely deter- 
mine the phase diagram of the system. 

Model. - In what follows (details will be given else- 
where [22]), we present a general formalism for a negatively 
charged fixed macroion interacting with Nc counterions, N+ 
positive salt ions, and N- negative salt ions in equilibrium 
with a bulk reservoir The Nc counterions are assumed to 
have a valency which may be greater than that of salt ions and 
therefore have to be treated differently. To proceed formally, 
we introduce the canonical partition function for the mobile 
charged species, confined to some arbitrary region around a 
fixed charged macroion characterized by charge density (T(r), 



Z = 



-H 



(2) 



Here, Vl{y) restricts the positions of mobile ions to the region 
of space available to them. The index a = {c, ±} stands for 
counterions (c), positive (+) as well as negative (— ) salt ions. In 
what follows, we assume that the macroion charge density a{v) 
is negative and confined to the macroion surface with surface 
charge density CTs. Introducing the density operator for each ion 
type Paijc) = ^{'^ ^ '^j)' Hamiltonian can be written 

in units of k-gT as 



W = 2 yy drdr'[gcPc(r) + g+P+(r) - q-P-{v) - cr(r)] 
X ^(r - r') [gcPc(r') + q+p+{r') - q-P-{v') - cr(r')] 
+ E[-y^"9X0)- / dr/i„(r)p„(r)], (3) 

where w(r) = ^B/|r| is the Coulomb interaction, and the gener- 
ating fields ha (r) have been added to calculate ion distributions 
by taking functional derivatives. Here we have also explicitly 
subtracted the infinite self-energies. 

At this stage, we proceed by applying the Hubbard- 
Stratonovich transformation [23], the purpose of which is to 
get rid of the quadratic density terms in H at the expense of 
introducing the fluctuating electrostatic potential field, (j){r). 
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This is followed by a Legendre transform to grand-canonical 
ensemble, where the number of ions {Nc, A^±}is replaced by 
their fugacities {Ac, Aj-}. Next follows the crucial step which 
makes the field-theoretic partition function convergent. We add 
the exponential of the following expression [21, 22] 



(4) 



to the partition function and subtract it perturbatively. Here 
UDH(r) is the inverse of the DH operator AttIq Vj^^{r — r') = 
(— V^ + K^) (5(r— r'), and of course corresponds to the screened 
DH interaction potential. Expressing the partition function in 
this form accomplishes three tasks: first, it removes the bulk 
density values of all ion types in order to make the one-particle 
densities finite. Second, the screening factor n^cfP makes the 
range of interaction between all the charges finite; and finally, 
the infinite self-energies cancel another set of divergencies 
present in the partition function. Rescaling all lengths by the 
GC length /.t = l/(27rZB'ZcO's) according to r = ^f, one ends 
up with an exact field-theoretic representation for the grand- 
canonical partition function Q~ J exp(— H/S), where H 
is the rescaled effective Hamiltonian [22] 

- ^q! [ drdf>(f)5DH(f-f')0(f') 
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2^ ^ dfO(f)g+(f)-— jAin{i)Q^{v). 

Here we introduce a shorthand Qi{Y) = exp[/i,;(f) — 
*<7,0(f) + (S/2)(g2/q2)i)oH(6)/2]-exp[(S/2)(g2/g2)A^;„] + 
(gf/2)(?i(f)2 - (S/2)(gf/(72)vj3H(6), while also rescaling the 
fugacities such that A^ = (27rS/i'^)Ac(. Here we also de- 
fined 47r[-V^ + (K/^)2]uj3jlj(f) 5{r - f'), and Auq = 
€'c(0) — wdh(O). The expectation values of different ion den- 
sities can be calculated by taking a functional derivative of the 
grand-canonical free energy with respect to the generating field 
hair), {pa{r)) = S\nQ/Sha{r)fi^\h^^n, giving rise to the 
rescaled densities 



(Pa(f)) = ^7^ = A„f2(f)(e- 
27r(Bcrf 



(6) 



where we have redefined A^ — *■ A^e ^ The nor- 
malization condition for the ion distributions then follows as 
Jdr[qcPc{r) + q+p+{r) - q-P-{r)] = qc J df CT(f), where 
a{r) = {p./ (Ts)<j{v). This corresponds to the overall elec- 
troneutrality of the system. 

Results. - Employing the Hamiltonian in eq. (5), we next 
make a full classification of possible limiting cases in terms 
of the coupling parameters {S, up); see fig. 1. These limiting 
cases are 



• (i) First, in the limit S <C 1, we find the familiar PB the- 
ory. This well-known regime [2,3,24] is characterized by 
many-body interactions among uncorrected ions. Math- 
ematically, it follows from the saddle-point equation for 
TtI [(/)], yielding the so-called Poisson-Boltzmann equation 
for the mean electrostatic potential 0pb, i-e- 



ql^^(t>PB 



(7) 



where = ±1 denotes the positive or negative sign 
of ions. Weakly correlated Gaussian fluctuations around 
the saddle-point solution may be captured by a loop- 
expansion in powers of S [17]. This regime separates into 
two sub-regimes according to the value of k/i, namely, the 
GC {kp 0) and the DH {up oo) regimes, which cor- 
respond to nonlinear and linear PB equations, respectively. 
The free energy of the system, F, is in both cases found to 
be a decreasing function of bulk salt concentration, but an 
increasing function of S, scaling as FpB ~ S [22]. These 
regimes have been throughly studied before [2, 3] and will 
not be considered here any further. 

• (ii) In the limit S ^ 1, we end up with the SC the- 
ory, representing a highly correlated system. Mathemat- 
ically, this follows from a virial expansion in powers of 
[17]. Since the electrostatic interaction between in- 
dividual mobile ions and the charged macroion dominates 
over the ion-ion interactions, the leading-order SC theory 
comprises a single-particle description with an effective 
ion-surface interaction potential of the form 



u{r) 



1 

2^ 



di'd{i') 



-K/j.|r — r' I 



-K^\v'\ 



(8) 

By adding enough bulk salt, the SC ionic distributions 
are destroyed, and the system ends up in the DH-regime, 
where electrostatic interactions are completely screened 
out. This shows up in the partition function as A± ^ 
{npY scaling [22], since the salt part of eventu- 
ally starts to dominate. The free energy is again a de- 
creasing function of salt concentration, but decreasing as 

The crossover from the PB to SC regime for the no-salt case 
has been extensively studied in the simulations [18], where the 
strong-coupling features are shown to set in at intermediate 
couplings about ~ 10. For the case with added salt such 
an analysis remains to be done. 

The asymmetric expansion of the partition function to the 
second order in Ac/S and to the first order in A±/E ^ 
is equivalent to the virial expansion used in the SC 
limit without added salt, together with the Mayer-Friedman 
resummation of the grand-canonical partition function for the 
simple salt ions, giving rise to the screened Debye-Hiickel 
potential [25,26]. Therefore we propose to call this expan- 
sion the Strong Coupling with Debye-Hiickel (SC-DH) theory, 
see fig. 1, identified already by Boroudjerdi et al. [1]. We 
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Fig. 1 : Schematic piiase-diagram representing different regimes of be- 
havior as a function of the two coupling parameters H and nfi. The 
crossover from PB to SC-DH regime takes place at sufficiently small 
salt concentration and by increasing the electrostatic coupling param- 
eter beyond a typical value of H, ^ 10 [18]. While, for elevated salt 
concentration nfj, > ~ 1, DH is the dominant regime. 



also expand different fugacities in powers of as Aq, = 
A° + Ajj/S + • • • , which allows us to avoid divergencies aris- 
ing from the second virial coefficient [17, 18]. 

Next we evaluate the ion densities to the lowest order in 
Ajj/S, and obtain ion density expansions as {pa) = + 
p^/S + • • • in the SC limit. This density expansion, listed be- 
low for all the ionic species (9)-(l 1), is in fact the main result 
of this paper. For counterion density we get 



eq. (4). These density differences are perfectly normalizable as 
demanded by the electroneutrality condition, whereas the den- 
sities themselves are not. 

The SC-DH can be evaluated in closed form only in the case 
of a single infinite plate, i.e. in the plane -parallel geometry 
with cr(f) = S{z) and counterions and salt ions being present 
on both sides of the plate. From eq. (9) we obtain to the leading 
order the SC counterion density as a function of the perpendic- 
ular distance z > from the wall, i.e. 

~pI{z) - Kl cxp { + [e--^~' - \]/np] (12) 

where the leading-order fugacity coefficient is given by 

A° ^np \ - 2{e(«++'?-)/'?="'^g+ + g_} + 

g,KMe9-/«=«^{Ei[g+/q,«M] - Ei[l/«:^] + log [q-/q+]} / 
2q_{Y.i[q+/q^Kp] - Ei[l/K^] + log [qc/q+W 
— 2q+e^''+^''^^/'''''^^ 

{-Ei[-g_/g,At/i] + Ei[l/K^] + log [q-/qc]} 



(13) I 



Thus, in the limit np counterion density approaches 
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iAi'n{Y')[e "('"^-e "(°°)] (9) This explicitly shows that we find the no-salt SC resuh in the 

limit fe- (/ ^ n n 71 i ^ ( 'A ^ Wr^t^a fVtcxf colt 
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and for salt-ion densities 



(p+(f)) =A°i^(f)e-^'^W+o(^[^|l! 



+ — u[r) 



o 



(10) 



(11) 



The second order terms can also be included in the above two 
equations, but they do not provide any further relevant in- 
sight [22]. Note that the key factor in the above expressions 
is the rescaled single-particle interaction term u{r), which cor- 
responds to a single ion interacting with the charged macroion, 
eq. (8). No assumption has been made thus far about the ge- 
ometry or symmetry properties of the macroion and the ex- 
pressions (9)-(ll) have a completely general validity in the 
SC limit. These formulas illustrate explicitly that we have an 
expansion in terms of the single-particle density differences 
|-g-ix(r) _ g-u(oo)j ^ [p(f) — p(oo)]. Stemming from the 
properly renormalized partition function with the counter-terms 



limit nfj, ^ [17], i.e. pl{z) — > e^^. Note that the salt correc- 
tion is negative for small distances, i.e., it reduces the density 
close to the charged wall. This is expected intuitively since 
the counterions tend to escape further away from the wall due 
to the reduced interaction in the presence of Debye screening. 
In the limit up oo, we find the usual DH expression, i.e. 

~ + c-'^^Vk/^), where p^l = Age-i/"'' is the bulk 
density. 

Therefore, the excess density effectively crosses over from 
one exponentially decaying form, i.e. the no-salt SC expression 
e^^//^ with the decay length, p (in actual units), to another ex- 
ponentially decaying form, i.e. the DH expression with the de- 
cay length K^^. In fig. 2 we show plots of the zeroth order ion- 
density /5°(z) relative to the contact value p°(0) for different 
values of np, clearly attesting to the fact that our results nicely 
interpolate between the well-known DH- and SC-regimes. 

The first order correction in 1/S for a planar charged surface 
follows again from eq. (9). The resulting expression can be 
given in terms of hypergeometric functions, but here we only 
present the result for np <^l, i.e. [22] 



z)- Kp{2i+l)\+0{[KpY), (15) 
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Fig. 2: Leading order counterion density in the SC limit for different 
values of nfi. For the sake of representation, we show the ratio of 
density to its contact value, pi{z)/ Pc (0), as obtained from eq. (12). 

which again exhibits a smooth transition to the no-salt case 
[17], and shows that also the first-order correction to the den- 
sity decreases close to the charged macroion surface when salt 
is introduced. Both, the zeroth- as well as the first-order counte- 
rion density profiles thus show a smooth transition to the no-salt 
case attesting to the consistency of our formulation. 

The concentrations of the salt ions can be obtained in the 
same way from eqs. (10) and (1 1) as 

= K\c^^{ + [e--^^^ -\]q+/{q,K^i)), (16) 

- K°_c^Y>{-[e-^^' -l]q-/{qc^^^i)], (17) 

for z > 0, which obey the overall electroneutraUty condition 
infinitely far away from the wall in the form 

This, together with the normalization, gives A*], ~ 

and Al = e^('J++'J-)/'J<:''^A']_, showing that negative ion den- 
sity vanishes exponentially fast as 0. It is not surprising 
that the positive salt-ion density shows a similar functional de- 
pendence on the distance from the charged wall as the counte- 
rion density. 

Let us now consider the validity of the SC-DH theory by 
comparing the leading order counterion distribution eq. (12) 
with the next leading contribution eq. (15). This gives a validity 
condition in terms of < 2S/(1 — K/i5), that makes sense as 
long as K/i < S^^. This means in fact that the SC-DH theory is 
valid for larger distances from the wall when compared to the 
zero-salt case [17]. This is clearly in accord with the features 
of the phase diagram in fig. 1, where by adding salt the SC-DH 
expansion eventually becomes valid for all separations from the 
wall. In the regime < k/z, this validity condition is always 
satisfied. In reality, this regime is easily applicable, since this 
condition can be written as q^nl-Q > 1. In physiological salt 
conditions with I-q ~ 0.7nm and k ~ Inm^^, this reads as 
0.7g^ > 1, which is already satisfied for divalent counterions. 

Finally, in the regime K/i > 1 we need to generalize the cal- 
culation to take into account corrections of order (k/x)^ to the 



next leading counterion density p\{z) in eq. (15). This calcula- 
tion will be presented in the forthcoming publication [22]. 

Furthermore, one has to notice that eq. (15) holds only in the 
case when the interaction between positive ions and negative 
salt ions does not give any significant contribution to the densi- 
ties. These plus-minus attractions induce an extra length-scale 
to the problem, the radius of ions a, which has to be non-zero 
to cut off these interactions. In the limit k,^ 0, i.e. when 
the interaction between ions is assumed to be unscreened, the 
results derive above hold as long as 

< ^ + 0[\n[ntB/{K.af]). (19) 

where a is the rescaled ion-radius [22]. This clearly means, 
that we cannot reduce a to zero without removing all the salt, 
i.e. setting also = 0. This is clearly caused by the fact, that 
even in the presence of very small amount of negative salt ions 
the Mayer functions of oppositely charged ions start to diverge, 
indicating complexation into Bjerrum pairs. 

In summary, the SC-DH limit, where typically one has a = 
a/ /J. 1, thus applies to the leading order for all values of z 
, as long as E!~^ < Kfj, < 2a/S. The criterion eq. (19) does 
not necessarily imply that the SCDH limit becomes invalid for 
larger values of K/^t. If exceeds the value given by eq. (19), 
one has to take into account the screening of interactions even 
between positive and negative salt ions themselves, thus mak- 
ing SCDH theory applicable for all values of salt concentra- 
tions [22]. 

The presence of these screened interactions increases the 
density of positive ions close to a charged wall and decreases 
the density of negative salt ions [22]. Intuitively, this is the case 
because single-particle density for positive ions is a decreas- 
ing and for negative ions an increasing function of distance 
from the wall. In the SC-limit, the plus- minus interaction is 
the dominant contribution to the second order virial expansion 
term. To maximize attraction between negative and positive 
ions all the density profiles become more sharply peaked close 
to the wall. However, this happens in such a way that elec- 
troneutrality holds, which furthermore means that the integral 
over the total charge density has to vanish to all higher orders. 
As a consequence, close to the wall, the amount of positive ions 
increases, and the amount of negative ions decreases. Clearly, 
the opposite happens far away from the wall. 

The second important contribution to next-leading order den- 
sities arises from the subtraction of artificial DH salt bath. This 
again gives positive contribution to positive ion density and 
negative contribution to negative ion density close to the wall, 
and is due to the fact that one is removing perturbatively the ar- 
tificial ions involved in the screening, thus decreasing the con- 
tact densities of negative ions and increasing the contact density 
of positive ions. 

The asymmetric virial expansion presented above also ap- 
plies to the high salt limit K/i ^ oo as long as S ^ 1. In this 
limit, the Mayer functions become extremely short-ranged due 
to an almost complete screening of electrostatic interactions. 
As a consequence, one can expand these functions in vdh, giv- 
ing rise to an expansion independent of S. This expansion is 
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valid only if WDH(2a) <C 1, where 2a is the ionic diameter. 
To the leading order, however, the density expressions derived 
above apply also to high salt case [22]. 

The present theory can be readily employed for further, more 
complex applications including charged polymers and colloids 
with biologically relevant concentrations of salt. One can also 
explore situations involving the joint interplay of counterions 
and salt within the present theoretical framework, which will 
be discussed elsewhere [22]. 

Another phenomenon that we have not yet accounted for thus 
far is overcharging and/or charge inversion observed in simu- 
lations of charged surfaces with multivalent counterions in salt 
solution [30-35]. This phenomenon has also been extensively 
studied in several recent theoretical works [16, 27-29, 36]. It 
is understood that ion-ion electrostatic correlations [16,29] as 
well as excluded-volume effects due to finite ion size [32, 36] 
play essential roles in the mechanism of overcharging. Physi- 
cally, the absence of this effect in our formalism so far may be 
traced back to the limiting one-particle structure of the leading- 
order SC-DH theory. A thorough analysis of higher-order terms 
is required to account for many-body and excluded-volume ef- 
fects which go beyond the scope of the present paper and will 
be presented in a forthcoming article [22]. It can be shown 
that the validity of the SC-DH theory strongly depends on the 
value of the ion-diameter through Mayer-functions, and can ac- 
tually break down in the limit where the interaction between 
oppositely charged particles exceeds k^T, which is already in- 
dicated by eq. (19). This limit indicates the onset of Bjerrum- 
pairing, thus destroying the Debye screening picture [22]. Our 
preliminary simulations [37] indicate that at large coupling pa- 
rameters overcharging is indeed possible but typically amounts 
to a fraction of the bare surface charge, which in this feature 
agrees with other simulations [31-35] and appears to reflect 
the second-order nature of this phenomenon. 

In summary, we have derived analytic expressions for 
counter- and co-ion distributions in the presence of salt in the 
vicinity of a charged macroion. The results are consistent with 
previous work in the PB, DH and SC limits and fill the gap be- 
tween the no-salt and high salt limits. The results presented 
here are relevant in a multitude of soft-matter and biologi- 
cal systems characterized by non-negligible salt concentrations, 
and pave the way for further applications [22], complementing 
those presented here. 
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